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Abstract 

The possibility of a magnetic phase transition in Heisenberg, Hub- 
K^ ■ bard, and s-f (Kondo-lattice) films is investigated. It is shown that, 

for any finite temperature {f3 < oo) and any finite number of layers 
^j- I (d < oo), the magnetization within every layer must vanish. Thus, 

^^ ' the Mermin- Wagner theorem is extended to a variety of system ge- 

^.j. • ometries. We also comment on the microscopic interpretation of the 

Q\ ', transition from two to three dimensions in the limit d,P -^ oo. 

1 Introduction 

C , In their well-known 1966 paper, Mermin and Wagner|l[] proved rigorously 

that, at any finite temperature, there can be no spontaneous magnetization 
in the one- and two-dimensional isotropic Heisenberg ferromagnet. They were 
also able to exclude the possibility of spontaneous sublattice magnetization in 
the one- and two-dimensional isotropic two-sublattice antiferromagnet. The 
proof is based on an inequality originally derived by Bogoliubov in a series 
of papers. [0] At roughly the same time as Mermin and Wagner, Hohenberg 
showed there could be no finite-temperature phase transition in one- and two- 
dimensional superfluid systems. His proof, too, is based on Bogoliubov's 
inequality. 
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Following these early papers, a lot of research has been carried out in or- 
der to extend the Mermin- Wagner theorem to other systems. Wegner (1967) 
considered a model describing a system with locally interacting itinerant 
electrons. 0] Walker and Ruijgrok (1968) discussed a band model for inter- 
acting electrons in a metal|P; Ghosh (1971), more specifically, recovered the 
Mermin- Wagner theorem for the Hubbard narrow-energy-band model. |0 A 
proof for the s-d interaction model is given by van den Bergh and Verto- 
gen (1974). A paper by Robaszkiewicz and Micnas (1976) extends the 
Mermin- Wagner result to a general model with localized and itinerant elec- 
trons, covering the modified Zener model, the extended Hubbard model, s-d 
models, and the model for the magnetic metal-insulator transition as par- 
ticular cases. 1^] In addition to extending the Mermin- Wagner theorem to 
several microscopic many-body models, generalizations to more complicated 
geometries were also found. E.g., Baryakhtar and Yablonskii (1975) prove the 
Mermin- Wagner theorem for systems with an arbitrary number of magnetic 
sublattices. Their proof not only excludes spontaneous sublattice magne- 
tization but also non-coUinear magnetic order with an external field being 
applied. Thorpe (1971) considers the case of ferromagnetism in phenomeno- 
logical models with double and higher-order exchange terms [p!0|. These re- 
sults were extended to the multi-sublattice case by Krzemihski (1976). [pJ| 
Recently, Matayoshi and Matayoshi (1997) have discussed models with n-th 
nearest neighbour exchange interactions and tried to extend the Mermin- 
Wagner theorem to anisotropic exchange interactions []T2|; their proof, how- 
ever, requires extremely special conditions on the parameters of the model 
which appear to be of no physical significance. 

The general procedure is essentially the same for most of the above papers: 
The Bogoliubov inequality is used with suitable operators defined in such 
a way as to give an upper bound for the desired order parameter, i.e. the 
(bulk or sublattice) magnetization. In general, one arrives at an upper bound 
that, in the limit Bq -^ 0, behaves as ~ {MBof/^ or ~ l/{\YiMBof/'^ thus 
excluding a finite value of the order parameter in one or two dimensions, resp. 
(M denotes the bulk magnetization and Bq is the external magnetic field). In 
three dimensions, no such behaviour is found, i.e. the usual line of reasoning 
fails and spontaneous magnetization cannot be ruled out. Within the usual 
scheme following Mermin and Wagner, the dimensionality of a system enters 
into the calculation only through the volume element when integrating the 
final version of the Bogoliubov inequality over fc-space. Thus, a microscopic 
picture of the transition from two to three dimensions does not emerge. 



Chester et al. (1969) argue|TB| that for three-dimensional Bose systems 



of finite cross section or thickness the results of the generically one or two- 
dimensional case are, in principle, reproduced. Their main line of argument 
rests on the operators in the Bogoliubov inequality being defined as Fourier 
transforms in a restricted space DxDxL or DxLxL, resp., where D 
remains finite in the thermodynamic limit while L goes to infinity. This way, 
the operators lack the physically intuitive meaning of the usual real-space or 
fc-space operators, which, again, means that no sound microscopic interpreta- 
tion can be assigned to the transition from two to three dimensions. Similar 
approaches have been used by Fernandez (1970) and by Fisher and Jasnow 
(1971) to discuss systems of distinguishable interacting quantum particles 
and spin systems, respectively. |T^ [|I^ 

Today, as a result of the rapid progress of thin- film technology in recent 
years, one can prepare and study systems with restricted geometries, such as 
films, in great detail. Thus, important parameters, such as the Curie tem- 
perature, magnetization and susceptibility, can be measured and discussed 
as functions of the number of layers d in a magnetic film. Experiments in- 
dicate that, for real systems, the transition from 2D to 3D behaviour of, for 
example, the critical exponent (3, occurs within a narrow crossover region of 
(j. |T^ The critical temperature Tc also shows a strong (i-dependence, with 



Tc{d) quickly approaching the bulk value Tc(oo) as d increases (see |]TB[,[0]). 

Taking these observations into account, it appears desirable to improve 
one's theoretical understanding of how the transition from the monolayer 
to the bulk occurs. By referring only to the truly microscopic properties 
of the respective many-body model, it should be possible to study in detail 
the transition from two to three dimensions. Work in this direction has 
been done by studying, both analytically and numerically, film systems with 
symmetry-breaking contributions to the Hamiltonian (for Heisenberg films 
with single-ion anisotropy see [@), thus allowing the study of the Curie 
temperature as a function of d and the anisotropy parameter. With respect to 
ferromagnetic itinerant-moment films, layer /overlayer geometries have been 
investigated |T^] , giving yet another indication of the relevance of many-body 
methods for film systems. 

At a more fundamental level, the validity of exact results, such as the 
Mermin- Wagner theorem in two dimensions, to thin films may be tested. In 
this paper, a proof of the Mermin- Wagner theorem for film systems within 
three main many-body models (Heisenberg, Hubbard, and s-f model) will be 
given. 



For reasons of simplicity, we shall restrict our attention to film systems 
composed of d identical layers stacked on top of each other. The calculations 
can easily be generalized to account for more complicated geometries. One 
can think of the film geometry as consisting of a two-dimensional Bravais 
lattice, i.e. the first layer, with A^ lattice sites and a rf-atom basis that 
corresponds to the d layers being stacked up. Every lattice vector then 
decomposes into 

Ria = Ri + fa (1) 

where Ri is a vector of the Bravais lattice and r^ is the basis vector pointing to 
the a-th layer. In all of the following calculations, Greek indices label layers 
and Roman indices refer to sites of the Bravais lattice. It should be noted 
that translational invariance can only be assumed within each layer. E.g., 
the notion of a reciprocal lattice only makes sense when referring to the two- 
dimensional Bravais lattice. A similar caveat applies to Fourier transforms 
being thought of as connecting real-space quantities with those defined in 
wave- vector space. 

2 Model Hamiltonians 

In the following we shall discuss the Mermin- Wagner theorem for the Heisen- 
berg model, the Hubbard model, and the s-f model. For film systems one 
has to distinguish between Bravais lattice indices and layer indices for all 
site-dependent quantities, such as spin operators S^^ -'^'V^^i ^ annihilation and 
creation operators c\J or coupling constants J°f which depend on two lattice 
sites. In this notation, the Hamiltonian for Heisenberg films is given by 

H = — 2_^ Jij {jia^j/3 + ^ia^jfs) ~'^Z_^^ '^ia (2) 

ijcxfi id 

where the term bJ2ia^~^^^'S^^ is due to the interaction with an external 
magnetic field b = ^■'^^ ° . The interaction with an external magnetic field 
leads to the magnetization 

where the phase factor e~^^'^^ accounts for both ferromagnetic and anti- 
ferromagnetic ordering, depending on the choice of K. We also assume the 



coupling constants Jj" to satisfy some general conditions, i.e. 

(4) 



ja/B _ jl3a 
ij ~ ji 


^/f = 




yRjn — Rpj 
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and 

/ ^ ^ X 2 

1 1 I Rm — Rn I 

= Q < OO. (5) 

■ye mp 

These conditions are very weak considering that for example the exchange 
integrals J will usually decay exponentially with distance. 

As an example for itinerant-electron systems, we shall discuss the Hub- 
bard model 



ijaj3a ia.u iot 

ija/3a ia ia 

where T^ describes the hopping of an electron from lattice site j of the /5-th 
layer to site i of the a-th layer, U is the energy associated with having two 
electrons at the same lattice site, and Ntot = Sjal^jat + ^jai)- Again, a term 
corresponding to an external magnetic field is included, crf^ in the Hubbard 
model is the z-component of the spin of the electrons associated with the 
respective lattice site, i.e. 



n 



0", 



ja 



(Cja|QaT ~ CjQ,|CiaJ (7) 



ial3 



Similar to the Heisenberg case, we require the hopping constants T'j^^ to sat- 
isfy the isotropy conditions T^^ = TJ^" as well as to converge upon summation 
over all lattice sites: 

,2 



In the s-f model one deals with two spin sub-systems, i 6*40, [ and {Bia}^ 
the former consisting of localized /-electrons, the latter being associated with 
itinerant s-electrons. With z t = ±1, the s-f Hamiltonian is 



H 



Erpa/S + 



ijaf3(7 



iP 



Za O; TLiacr + >^icyCi„_„Ciacr 



i+^e-^-^»(5^+a^ 



(9) 



3 Bogoliubov inequality and choice of opera- 
tors 



The original Bogoliubov inequality as derived in P[ is 



\C,Al 



- 2 



A,A^ 



[C^HUC 



(10) 



where A and C are local operators and (....) denotes the thermodynamic 
expectation value. It should be noted that the two factors on the r.h.s. each 
are, mathematically, upper bounds to a norm and, thus, greater than or equal 
to zero. In particular, if, for example, the double commutator depends on 
some parameter fc, we will always have 



[C,HUC+] ){k)+{\[C,HUC^] ){k')>{\[C,HUC^] ){k){ll) 



For our purposes a slightly modified version of the Bogoliubov inequality 
shall be used. Dividing both sides by the double commutator and summing 
over all wave vectors k associated with the two-dimensional Bravais lattice, 
one arrives at 



[C,Al 



?([[C,i7]_,C+]_ 



^fE([A^ 



(12) 



The choice of suitable operators A and C is crucial; it determines whether 
the inequality will be physically meaningful or not. In film systems long-range 
magnetic order within a given layer is conceivable where the bulk (or even 
sub-lattice) magnetization of the whole system vanishes, e.g. 

• •• Tinnnnnnn... 

• •• nnnnnnnn--- 
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Excluding long-range magnetic order for every layer within a film would, 
therefore, be a considerably stronger statement. The general idea of our 
proof is to use the Bogoliubov inequality to find an upper bound for the 
layer magnetization Mp, i.e. 

M^</(i?o,M) 

where / is a function that approaches zero as i?o ^ and does not depend 
on any layer-specific quantities. 

This is best achieved by choosing C as a sum I]q,(...) of spin operators 
and A as a spin operator associated with a specific layer (3, say. This way, the 
numerator of the l.h.s of the Bogoliubov inequality will be a layer- dependent 
quantity, while the double commutator in the denominator will be summed 
over all lattice sites and thus be layer-independent. One may then expect to 
be able to replace the r.h.s. of the inequality by a (layer-independent) upper 
bound. More specifically, for the Heisenberg model we set 

A(^^) = S^{-k~K) (13) 

and 

C^ECp^ESMk) (14) 

where we have used the Fourier transform 

S^+'-'^'y''\k) = J2 e^'-'^^Si^'^'"'"''^ (15) 



Note that in /c-space we have (S'j^{k)] = Sp{—k) due to the definition of the 
Fourier transform. The operators S in the context of the Heisenberg model 
are spin operators for which the usual commutation relations hold true, such 
as, in A;-space, 

'Si{h),S^ih)]_ = 2h5^f,S:,{h + k2) (16) 



and 



S'^{h),S^{h)]_ = ±n6^f,S^{h + k2) (17) 



In the Hubbard model a similar definition for the operators A and C is 
used. However, as mentioned above, the operators S**^"-* now are built up 
from fermionic creation and annihilation operators, e.g. 

^.-(^i) = H.i^^T (18) 



and 



(19) 



(20) 



(21) 



respectively. The commutation relations for spin operators may, in a purely 
formal sense, be used in this case as well. 

In the s-f model, one must be careful not to forget that one is dealing with 
two separate spin sub-systems, one of which can be described by usual spin 
operators, while the other is associated with itinerant electrons. The two 
systems are independent from one another in the sense that spin operators 
and creation or annihilation operators commute: 









a+{k,) = ncl^^c^,^^^ 






The operators 


A and C are, then. 












A(a) 


= ( 


y-A-k- 


K) = 


-k-K,c 


A^-k- 


KM 


and 










j 




\ 




C = 


/3 


a^ik) = 




\ 


= >: 

/3 


cA 



s. 



{+,-,x,y,z) (+) 







and, thus, 



s^^'-'-y-\h),c^^X. 



(22) 
(23) 



With (|TP|),(|T^) we define as operators A and C in the s-f model 



A 



(7) 



S-{-k-K) + a-i-k-K) 



and 



/3 /3 



(24) 
(25) 



4 Evaluation of the Bogoliubov inequality 



4.1 Hamiltonian- independent quantities 



It is obvious from the structure of the Bogoliubov inequality ([TJ) that the 
numerators on both sides of the inequality are determined entirely by the 
choice of operators as discussed in the previous section. The details of 



the many-body model enter the calculation only via the double commutator 
[[C, H]_,C^]) . For simplicity, we shall, therefore, start with the quantity 



[C,A]_ 



J2/[s;{k) + a;{k),S-i-k-K) + a:;{-k-K) 



(26) 
(27) 



mnf3 



Jk-Rm^-i{k+K)-R„ 



^mf3 + '^'^mfSl^rnPi, ^n-y + '^'^n-yi'^nfl 



for the s-f model. For reasons described above, S and a operators com- 
mute, so the commutator can be evaluated directly using the fundamental 
commutation relations. This leads to 



[C,A]_ 



mn/3 



where we have introduced the layer magnetization 

M,{T, 5o) = ^^ E e^^^"^" (SL + < 



N h 



n7 ' ?17 



(28) 



(29) 



It is obvious that in both the Heisenberg and the Hubbard model the 
same relation 



C,A 



(7) 



-M^{T,Bo 



(30) 



gjfj'B 

holds, the only difference being the fact that only one spin system contributes 
to the magnetization. 

The r.h.s. of the Bogoliubov inequality is proportional to the anticom- 
mutator sum 

E(Ka+]^) (31) 

k 

which, for the s-f model, is given by 



E[AA^ 



kmn 



'i(k + K)iRm-R-a) 



^n-f ' "''^n7j,'^"7T 



'^m,7 ' '^^■m-i'\^rn-yi 



(32) 



Summing the exponential over all k gives the delta function NSmn- The mixed 
commutators involving both spin and creation/annihilation operators in this 
case do not vanish, since we have the anticommutator rather than the com- 
mutator. Thus, we arrive at the expression 



E[AA^ 



NT. 



O- 0+ 



+ Zll ( C^^iCn-yfO^^ + C„^|Cn^jD„^ 



"T"" \'^n7i'^n7T'^n7t'^'^7i "■" '^n^^'^n^i'^n-yi'^ri-y^ / j 



(33) 



Since we are interested in the quantity (J2k[^^ ^^] + ) ^^ ^^ upper bound 
in the Bogoliubov inequality, it suffices to find upper bounds for the individ- 
ual expectation values on the r.h.s. For the Heisenberg term we find 



/ . \ |_*-'n7' ^n^ 



2 E {iS:,r + iS'n.r} < 2 E {si,} = ^^'S{S + 1)N 

n n 

(34) 
For the Hubbard contribution we have 

n n 

< E((^«7i) + (^«7T)) 

n 

< 2N (35) 

For the mixed terms appearing only in the s-f model, an upper bound is 
given by 

(36) 



T «,c„7-^^n7> < 2 (4 + 2SiS + 1)) N 



as shown in appendix A (see eqn. (|7TD). 
Tabulating these results we have 



^/[a,A+] \ <2h^S{S+l)N^ 
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for the Heisenberg model; 



e(k^^ 



< 2h'^N^ 



for the Hubbard model; and 



5]/[A^+] )<N^h^ASiS+l) + lO) 



for the s-f model. 



4.2 The double commutator ([[C, //]_, C+] 



In this section, we shall calculate the remaining double commutator ([[C, H]^, C^ 
for the individual models and give upper bounds for use in the Bogoliubov 
inequality. 



4.2.1 The Heisenberg case 

The double commutator ([[C, H]^,C^]) in this case is 



I 7e \ '- 



■ye mp 



,-ik-{Rm-Rp) 



,S;i-k) 



(37) 



The real-space commutator on the r.h.s. can be easily evaluated using the 
standard commutation relations. Thus, one arrives at 



[c,HUc^])ik) = T.T.Jpl^Hi^-^"'-^'^-''^^^ 



7e mp 

■ /2S'^ S*^ -I- 9+ 9^ 
\ perri'y '^ ni'y pe 



(38) 



-2bh'j:j:e 



-iKR„ 



Krnel 



To this we add the double commutator ^[[C, iJ]_, C+]_ ) (—A;) which, as dis- 
cussed above, is a positive real number. Replacing the spin operator expec- 
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tation values by the upper bound 2h S , we find, after some minor algebra, 



7£ rap 



[C,HUC+\_j<ANdh'\BoM{T,Bo)\+8h^S'J2T.\jpl 
where we have already used the fact that 



1^ ■ [-tim J^p 



4 



(39) 



1 — cos(A; ■ {Rm — Rp)) < 



1^ ' \Rm Rp 

4 



(40) 



With the above definition of the constant Q, we arrive at the final result 



[C,if]_,C+]_\ < ANdh^ {\BoM{T,Bo)\ + 2h^S^Qt 



(41) 



4.2.2 The Hubbard case 



Again, we need to calculate the full double commutator ([[C,H]_,C~^] 
With the Hubbard Hamiltonian given in eqn. (BI), we have 



T.Hik),H 



(42) 
which still needs to be commuted with Y,s^e{~^) ■ Replacing the expecta- 
tion values by their modulus, one gets the relevant inequality 

[C,HUC+]){k) = J2(\[^^ik)^Ho + H,]_,a;{-k) 



< h^Y. ^r (e-'^-(^'-^') - 1 



ilae 



+2bh'j:(<iK) 



CjotQeT 



For the same reasons as above, we may now add ([[C, i/]_, C"*"] 
get an upper bound for the l.h.s. For the expectation values 



'^lei^ial 



(43) 



-k) to 



Ct^'^t 



+ 



C;_ I C. 



lei'-iai 



Cicc^Cle] 



Ci^ I c. 



lei'-iai 



< 4 (see 



we find as an upper bound 
appendix A, eqn. (|7^ , for detailed calculation), so in total we may write, 
using the above notation. 



[C, H]_, C+] ) < ANdh^ {\BoM{T, Bo)\ + 2qk^) . 



(44) 
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4.2.3 The s-f case 

The s-f double commutator 

' 7e 



S^{k) + a+{k), Ho + Hh , S;{-k) + a;{-k) 



(45) 

may be computed by considering first the field-independent contribution to 
the Hamiltonian 



^0 = ETi 



a/B 



J 



ija/3a 

and the respective commutators 



ij ^iacr^jl3cr ^ 2_^ (^^o-'-'ia'^icwT + '^ia'^ia-a'^ic 



(46) 



^ml3^^0 _ - ^"^krrfi 



^mlS'i^rnPl '^k-yl'^mf3l 



jn( Strata -2SL^a, 



mfi'^rnP '^ '-' mP'^ mfi 



k-y 



'S'm/3) Hq 



J ri b^aCT^o Zb, 






(47) 
(48) 



We then find (neglecting, for the moment, the external contribution Hj)) 



kmn P^S 



'^mlSl'^nifSl C^7tCm/3J, I , I O^^ + C^Sl'^nS'l 



[ 

kn (3^5 



(49) 



Following the usual procedure of adding ([[C,H]^,C^]) {—k) we arrive at 
the upper bound 



[C,Ho].,C' 



nA;7/3 



1 — COS (k-(R„~R 



< 2h'j:\T:^\k'{Rn-Rk 

nfe7/3 

= 2Ndqn^k'^ 



(50) 
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where (c^/3|Cfc^j,) + (c^/jiCfc^t) ^ 4 has aheady been used. 

To this prehminary resuh, we still need to add the double commutator 
([[C, Hh\-,C^]_) with the external part Hi, of the Hamiltonian. Here, too, 
we need to take into account the presence of two distinct spin systems. Thus, 
Hf, is given by 

H, = -Bo^ Y.i9jSL + 2arJe-^-^» (51) 

ia 

where the usual definitions of the spin operators apply. As in the previous 
cases, one gets 



[C,H,UC' 



-2h^NdBoM{T, Bo) 



where, however, the magnetization is now defined as 

1 fJ'B 



M(T,5o 



Nd h 



Ee-^^-^»((^.5f,)+2(af, 



i(3 



Finally, we arrive at 

'[C, //]_, C+] ) < 2Ndn^ (\BoM{T, Bo)\ + qk'^ 



(52) 



(53) 



(54) 



5 Proving the absence of spontaneous mag- 
netization 



Within all three models, we have found for the double commutator 



[C, H]_, C+] ) < eoNd (\BoMiT, Bo)\ + iik 



(55) 



where ^j are constants depending, at most, on fixed parameters of the re- 
spective many-body models. We also know that 



and 



[C,A]) = i^NM,{T,Bo) 



E(K^1+)<2e3iv 



(56) 



(57) 



We can now give a generally applicable discussion of the layer magneti- 
zation. With the Bogoliubov inequality (|I2D, 
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[C.Al 



2 



/9 



V l\ ' , <^V(k^+1 ) (58) 



we find 

Y gM (|BoA^(r, B„)| + fty) -^ 

The l.h.s. of the inequahty can be replaced by an integral, using the formula 

where j^^ is the area in two-dimensional /c-space associated with one quan- 
tum state. Restricting the support of the integral to a finite-volume sphere 
inscribed into the first Brillouin zone only strengthens the inequality, so 

where /cq is the cutoff corresponding to the sphere in A;-space. In the ther- 
modynamic limit, "^ = "^0 approaches a constant finite value as A^ ^ oo. 
Evaluating the integral and performing some minor algebra, we have 

M'JT, Bo) < ^— ^^,,. , (62) 




' In (l + 



gifcg 



\BoM{T,Bo)\ 



(^ again is an unsignificant constant). 

From this formula, it is obvious that, for any finite temperature {(3 < cxd) 
and finite thickness {d < oo), the logarithm in the denominator will diverge 
in the limit Bq ^ 0, thus forcing the layer magnetization M^ to vanish. 
One should note that the final result does not depend on the choice of K 
and, thus, excludes both ferromagnetic and antiferromagnetic ordering. This 
proves the Mermin- Wagner theorem for film systems within the Heisenberg, 
Hubbard, and the s-f (Kondo-lattice) model. 
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6 Summary 

In this paper, we have been able to prove that the Mermin- Wagner theo- 
rem which was originally shown to exclude a magnetic phase transition at 
finite temperature for one and two-dimensional systems can be extended to 
Heisenberg, Hubbard, and s-f (Kondo-lattice) films of any finite thickness. 
While this may be expected from very general considerations, and similar 
theorems for Bose systems and essentially classical systems were proposed 
before, a microscopic calculation for many-fermion systems has been lacking 
so far. 

In view of recent progress in the experimental and theoretical study of 
thin magnetic films, it appears to be worthwhile to reconsider the impor- 
tance of the Mermin- Wagner theorem. This is particularly the case when the 
symmetry-breaking is not due to explicitly adding anisotropic contributions 
to the Hamiltonian, but is instead a secondary result of using numerical ap- 
proximations that violate the Mermin- Wagner theorem. Future work in this 
direction might be able to provide useful insights for theorists and experi- 
mentalists alike. 

Our microscopic approach suggests that the parameter d, i.e. the num- 
ber of layers in a film, plays a similar role as the inverse temperature j3. If 
either one diverges, a phase transition cannot be ruled out. It is, of course, 
not a priori clear if the behaviour of an upper bound on a physical quan- 
tity has any physical significance itself; however, we believe our results to 
indicate that, within the discussed models, in order to possibly describe 
(anti)ferromagnetism, one has to take the thermodynamic limit seriously, 
in the sense that N ^ oo does not suffice, but rather c? ^ oo is required as 
well. 
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A Upper bound for E^ (c^ Cn^-aS^^) and 



c 



+ 



n 



.fja^k-fa 



We first discuss the quantity J2a \'^i'ya'^n'y-crS~^). The problem in evaluating 
the expectation value arises from the fact, that we are dealing with seemingly 
unrelated operators. It is, therefore, convenient to use the identity 






~'^ y^nycr'^ny-cr ~ '^^nyj \^ny-a^nya + ^'-'717 j J 

= \ll<^ij)B,<rB% (63) 

where j labels the individual terms and is a phase factor (0(1) = +1, 
0(2) = —1, 0(3) = i, 0(4) = —i). The original expectation value has thus 
been decomposed into the sum of expectation values of pairs of adjunct 
operators Bj„B'^„. Using the spectral theorem, we find 

E {<.Cny^.S-;) = ^ E / dE^^^^^m E S^^,^^ (E) (64) 



where S + (E) = ^ ( -B^, Bj^ ) (E) is the spectral density in its energy 

representation. The l.h.s. is the expectation value of the sum of two adjunct 
operators and, thus, is real. We may now make use of the fact that for pairs 
of adjunct operators, the spectral density is positive definite, which, together 
with the triangle inequality and |0(j)| = 1 gives an upper bound for the r.h.s. 



(65) 
The sum on the r.h.s. now consists of the 0-th spectral moments associated 
with the operator pairs Bj'^, Bj^. Each of the spectral moments is given by 



the expectation value of the anticommutator 

r(0) 



^ja^3<y 



— oo 

The anticommutators can be easily evaluated: 



Bt,Bja 



(66) 



M 



(0) 






M 



(0) 

Bo„B2a 



< 4 + 25(6 + 1) 2 ( D c^ ^^c„^cr + 6 c c. 



^n7 '^nya'^ny-a 



M 



M 



and hence 



(0) 
(0) 



< 4 + Zb\b + I) ~\- Zl (bC^^Cn-ycr ^n-y ^n-icr^ny-a 



< 4 + 2S{S + 1) - 2^ Isic+c^,, - S-; 



^n7 '^n7o-'^n-7-o" 



E (c+^,c„,_.5„7) < 2 (4 + 25(5 + 1)) 



(67) 
(68) 
(69) 
(70) 

(71) 



For the expectation value 



one uses the same procedure. For 



{n, P) = (A;, 7) it is obvious that \Cni3a^n(ia) < 1- Starting with the decompo- 
sition 



^iaa^li"^ 



~ [^iaa ~ ^Iva) {'^iaa ~ Qua) 
+^ {'^iaa + '^'^luaj {Qaa '~ IClva) 



^ y^iaa '^QuaJ [Qaa + l-Cluajj 



(72) 



one can write the expectation value as 

4 
4 



^nfia^kya 



1 4 1 °? -, 

3 —00 



— A '^ AiA 



i=l 



J J 



(73) 
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(similar notation as above). 

Since we have discussed the case (n, P) = {k, 7) separately, we can now 
assume (n, /9) 7^ {k,'^). The 0-th spectral moments can be easily calculated; 



they each give M 



(0) 



2. Thus for all n, A;, /3, 7 we have as an upper bound 



^nlBa^k-ya 



< 2. 



(74) 
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